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INTRODUCTION 


Ultrasonic  measurements  at  very  high  pressures  give  us  the 

unique  opportunity  of  obtaining  the  elastic  properties  of  crystalline 

solids  experiencing  extremely  large  elastic  deformations,  which  cannot 

be  attained  by  other  more  conventional  techniques.  At  best  conventional 

stress-strain  measurements  in  crystalline  solids  indicate  elastic 

strains  of  the  order  of  one  percent;  in  sodium  chloride  the  linear 

strain  at  a hydrostatic  pressure  of  270  kbar  is  nearly  eleven  percent. 

Methods  of  continuum  mechanics  brought  to  bear  on  high  pressure 

problems  were  inspired  by  Bridgeman's  pioneering  experimental  work. 

2 

Murnaghan  did  the  early  theoretical  work,  inventing  the  finite  strain 
method  and  using  it  to  third  order  elastic  constants  for  the  equation 

3 

of  state.  Birch  was  another  early  worker  using  a somewhat  different 
though  valid  formalism.  A treatment  giving  the  fourth  order  elastic 
constants  for  single  crystals  using  lattice  dynamics  was  given  by 

4 

Ghate  . The  present  treatment  derives  both  the  equation  of  state  and 
the  velocity  in  terms  of  the  fourth  order  elastic  constants  by  methods 
of  continuum  mechanics. 

The  methods  used  here  are  based  on  the  tensor  formalism  developed 
in  Green  and  Zema.^ 

^Green,  A.  E.,  W.  Zema,  Theoretical  Elasticity,  Oxford  at  the 
Clarendon  Press,  1954. 

^Murnaghan,  F.  D.,  Finite  Deformations  of  an  Elastic  Solid,  Dover 
Publications,  Inc.,  New  York. 

^Birch,  F.,  "The  Effect  of  Pressure  Upon  the  Elastic  Parameters  of  Isotropic 
Solids,  According  to  Mumaghan's  Theory  of  Finite  Strain,"  Jour.  Appl. 
Phys.  9,  279,  1938. 

^Ghate,  P.  B.,  "Fourth  Order  Elastic  Coefficients,"  Jour.  Appl.  Phys. 

35.2,  337,  1964. 


The  experimental  techniques  for  obtaining  high  "quasi-hydrostatic" 
pressures^  and  for  making  high  pressure  ultrasonic  interferometry 
measurements^  were  described  in  detail  elsewhere  and  are  shown 
schematically  in  Figure  1. 

The  direct  experimental  data  obtained  from  ultrasonic  interfer- 
ometry consists  of  the  frequency  intervals  Af_  and  Af_  (for  dilata- 

P * 

tional  and  shear  modes  respectively)  for  maximum  destructive  interfer- 
ence of  the  first  echo  of  the  ultrasonic  echo  train  as  observed  on  the 
cathode  ray  tube.  Here  the  respective  velocities  are  given  by  (Af)/(2d) 
where  d is  the  specimen  thickness.  The  isotropic  velocities  were 
obtained  through  the  use  of  the  Decker  equation  and  the  frequency  inter- 
vals in  the  universal  identity  relating  the  adiabatic  density  deriva- 
tive of  the  pressure  and  the  dilatational  and  shear  isotropic  velocities 
Vp  and  Vg  respectively  at  all  hydrostatic  pressures:  ^/^Pad  “ ^p  * " 
(4/3)Vs^ 

DERIVATION  OF  VELOCITIES  FOR  THE  THEORY  OF  SMALL  DEFORMATION  SUPERPOSED 
ON  FINITE  UNIFORM  DEFORMATION 

Using  the  notation  of  Green  and  Zema^  and  conforming  to  the  exper- 
imental situation  we  apply  a uniform  hydrostatic  deformation  to  the 
unstrained  body  Bq  which  becomes  B in  the  new  (pressurized)  state.  To 
the  body  B we  apply  an  infinitessimal  deformation  (the  ultrasonic  wave) 
at  time  t,  and  represent  this  configuration  as  B'  (Figure  2). 

^Green,  A.  E.,  W.  Zema,  Theoretical  Elasticity,  Oxford  at  the 
Clarendon  Press,  1954. 

^Kendall,  D.  P.,  P.  V.  Dembowski,  T.  E.  Davidson,  Rev.  Sci.  Inst. 

Vol.  46,  No.  5 (1975). 

^Frankel,  J.,  J.  F.  Rich,  C.  G.  Homan,  J.  Geoph.  Res.  Vol.  81,  No.  35, 
6357-6363,  1976. 
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Taking  the  moving  coordinates  to  coincide  with  a fixed  rectangular 
system  {x,y,z)  in  the  strained  body  B (Eulerian  formulation)  we  have 
©1  = X ©2  = y ©3  = z (strained  body)  (1) 

The  rectangular  axes  which  define  points  in  Bq  are  taken  to 
coincide  with  the  axes  (X2,X2,X3). 

The  constant  extension  ratio  Xj,  (the  strain  parameter)  is  given 


by 


Xi  = JL=  JL=  ^ 

XI  yj  zj 


(2) 


By  taking  the  body  (working)  coordinates  in  the  deformed  body 
(l.e.  the  working  coordinates  move  with  the  mass  points  of  the  body) 
and  using  the  definition  of  the  metric  gj^j  from  the  equation 
ds*  = gjj  d©^  d©^  we  have  the  metric  tensor  of  the  undeformed  body 


and  g = l/Xj® 


Here  the  sub  or  superscripts  refer  to  the  covariant  or  contra- 
variant  tensor  representations,  g is  the  determinant  of  gj^j.  In  the 
same  representation  the  metric  tensor  of  the  deformed  body  B is 

= G^j  = ; G = 1 (4) 

where  6 is  the  Kronecker  delta. 
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The  Deformed  State  (Under  Finite  Deformation] 


The  strain  invariants  are 


= 3X, 


U = 3X, 


I3  = X,‘ 


and  the  strain  tensor  elements  become 

^ij  ~ 1 T ■ ^^2)  (6) 

The  strain  invariants  (Eq.  5)  are  independent  of  the  coordinates, 
they  are  constant  for  a given  deformation.  At  this  point  we  invoke 
Murnaghan's  theorem  (p.  61  of  reference  2)  for  isotropic  elastic  media: 
"A  deformable  medium  is  isotropic  if,  and  only  if,  the  strain  energy 
density  is  a function  of  the  three  strain  invariants  Ij,  I2,  I3". 

We  can  now  express  the  strain  energy  density  W in  terms  of  Ij, 

I^,  and  after  some  manipulation  (as  given  in  reference  2)  and 
using  the  principle  of  virtual  work  we  can  obtain  the  components  of 
the  stress  tensor  of  the  body  under  finite  deformation. 

= ^Xj*  + 24'Xj'*  + p 

= 0 (71 


where 


$ - M 


y _ 2 aw 
■ ^312 


= 2/T7 


It  should  be  noted  here  that  since  the  unit  reference 

area  is  taken  in  the  deformed  coordinates.  Equations  (7)  therefore 
represent  the  physical  components  of  the  stresses  and  in  particular 
T^^is  the  hydrostatic  pressure  in  our  experiment  (t^^=  -P), 


^Mumaghan,  F.  D.,  Finite  Deformations  of  an  Elastic  Solid,  Dover 
Publications  Inc.,  New  York. 
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The  Perturbed  State 


We  are  now  in  a position  to  derive  similar  quantities  for  the 
perturbed  case  B' , the  displacements  due  to  the  ultrasonic  wave. 

In  this  new  configuration  the  displacements  can  be  written  as  the 
sum  of  the  initial  uniform  displacement  of  the  body  v and  that  of 
the  perturbation  ew,  v (0i.02»®3» ^w(0i,02»03.t)  where  e is  a 
constant  so  small  that  the  squares  of  e and  its  higher  powers  may  be 
neglected  compared  to  e. 

The  covariant  base  vectors  of  the  coordinate  system  0^^  at  points 
P'  of  B'  are  denoted  by  Gj  + eGj  = r 4 + v i+ew^  where  r is  the 

J.  X ,X 

location  of  point  Pq,  and  the  commas  indicate  differentiation  with 
respect  to  the  i^^  coordinate. 

To  the  first  order  in  e the  perturbed  metric  tensor  C^he 

perturbed  state  is  usually  denoted  by  primed  quantities)  of  the 
deformed  body  referred  to  base  vectors  of  body  B is: 


"ij  = 


j 3u 

3u  ^ 3v 

3u  3w 

3x 

3y  3x 

X 

ro 

CO 

3v  ^ 3w 

(symm) 

2 — 

3y 

3z  3y 

2 ^ 
3z 


- G 


ij' 


(9) 


and  ojj  = u 012  = v 


0)3 


= w 
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The  perturbation  in  the  strain  invariants  is 


1 1 3y  3z 

I-,  = 4X,“  (in.  + li  + 

^ ^ 3x  3y  3z 

I . = 2X  ® f—  + — + 

^ 3 ^^1  ^Sx  ^ 3v  3z5 


(10) 


^3x  3y  3z‘^ 

The  perturbation  components  of  the  stress  tensor  are  given  by 
(see  ref.  [1]  p.  119) . 

,11  _ P 9u  p 

^ ^11  37  ^12  37^ 


t.22  , C ^+Co(^+^) 
^ 11  3y  3x  3Z'' 

t-33  = C,,  Ci2  ^ 

^1  3z  3x  3y 


(11) 


t' 

|12  , 

C44(|H 

+ IV) 

3x 

t' 

,23  ^ 

^44  (|l 

+ ^ 

3y 

t' 

,31  ^ 

^44 

+ ^) 
dz^ 

Where 

^11  = 

CM 

+ 

' + 8BXj® 

12  10 
+ 2CXj  + 8DX1  + 

II 

CM 

u 

-Xi^, 

j)  + p + 

ZAXi**  + 

8BX1*  * 2^^1^ 

-12) 


(12) 

1®  + 8FXj® 

i[°  + 4EXj®  + 8FX1® 
(13) 


J_i!w  ^ B =-i-i!w 
/I7  3Ii*  ’ /ij  312"  ' 


Q _ _L_ 

■ /I7  313^ 


D = 


3^W 
/TJ  3I23I3 


E = 


3*W 


/TJ  3I33I1 


F = 


2 3^W 

A7  3Ij3l2 


(14) 


^Green,  A.  E.,  W.  Zema,  Theoretical  Elasticity.  Oxford  at  the 
Clarendon  Press,  1954. 
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r s 

The  physical  components  of  perturbation  in  stresses,  t'  , can 
be  obtained  by  the  contravariant  transformation  rule,  and  to  the 
first  order  in  e and  referring  to  the  B'  configuration  we  have 


^'rs  ^ ^’rs  + ^ms  ^ ^^rm^ 

80™  90"' 


(15) 


Using  (12)  we  have 
'll 


, |1 1 , o-rll^  9u  . p /-Bv  . 9wn 

t’  - CC„  * 2t  ) * Ci2(—  . _) 


^(22  _ fft  ^ 11\  3v  . p ^3u  . 9 


t.33=  (Cii  . 2x^1)|^.C,2(|H.|^) 
t.23=  (C44.t11)(|I>|^) 


(16) 


t.^l=  (C44.xlb(|^.|H) 


The  equation  of  motion,  with  f as  acceleration  components,  become 


3t' 


30‘ 


= pf’^ 


(17) 


From  (17)  and  (16)  we  have  (after  some  manipulation) 

il. ,3*u  3*v  3*w 


+ ^ + — ) = P 772 


3x^  3x^  3x3z" 


3^u 

3t= 


T (Cii-Ci2+2x^^)V2v  + j (Cii+C,2+2x^1)(^  + * ^^)  = p 

2 ^1  12  ^ 11  12  3y2  3y92  3x3y 


3fv 

3t* 


i (C„-C,,-.2x^^V^w  i (c„+C,,-.2x^^)(^+  ^)  = p ^ 


2 ^"11-^12^ 


V'-l  ,"‘•1  J -r  J 

2 12  3z2  3x3z  3y3z 


3t" 


(18) 
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Vectorial ly  adding  above  with  w = we  have 

J * 2 w)  = p w (19) 

This  is  the  equation  of  motion  of  body  B'  in  terms  of  perturbation 
(w)  of  displacement  of  body  B. 

Taking  the  divergence  and  curl  of  (19)  we  have 

THE  UNIVERSAL  IDENTITY 
From  (16)  we  have 

^,11  + t'22  + t'”  = (Cjj+2Cj2*2t^^)  (I7  + I7  + It)  (21) 


'■3X  3y  3zJ 


t'll  = t'^^  = = - A (Pressure)  = - AP 


9u  ^ ^ _ - A (Volume)  _ ^ 

3x  3y  3z  ~ «.  V j ~ ^ 

Unperturbed 

Volume 

hence  ^ = - A?  (23 

P 

(23)  and  (22)  with  (21)  give 

p — _ Limit  p — = i (Cl  1+2C,  ,+2t^^)  (24 

dp  AC  0 Ap  3 

Hence  from  (24)  and  (20)  eliminating  have 

V * - — V * + (25 

P ' 3 s dp 

(25)  is  the  universal  identity^* 

^Truesdell,  C.,  Handbuch  der  Physik,  EdLS.  Flugge,  Vol  III/3,  Springer 
Verlaft  Heidelberg,  New  York. 


From  (20)  we  have  the  velocities 


Vp=  C- 


Cn+2T^^ 


Dilatational  Velocity 


Cii-Ci2+2t^^  1/2 

V_  = (—Li — — Shear  Velocity 

2p 

Where  C^j,  etc.  can  be  computed  in  terms  of  derivatives  of  W with 
respect  to  strain  invariants  from  (13),  (14)  and  (8). 

APPLICATION  WITH  FOURTH  ORDER  ELASTIC  CONSTANTS 
The  Velocities 

In  this  section  we  assume  that  the  elastic  strain  energy  for  the 
material  under  consideration  can  be  expanded  up  to  fourth  order  terms 
in  strain  invariants.  The  strain  invariants  defined  by  Mumaghan^ 
are  slightly  different  than  those  of  Green  and  Zema^.  For  the  sake 
of  easy  comparison  we  give  the  transformation,  denoting  those  of 
Murnaghan  by  a superscript  c, 

i (Il-3) 

^2  " T •^^^2-^5  -2(Ii-3)]  (27 

I3  = i [(I3-I)  ♦ (Ii-3)  - (I2-3)] 

W = j(X+2y)  if  - 2y  I^  + jC*  . 2mljl=  n I^ 

+ 16q(l5)'‘  + 16r  (ij*)  (I^)  -Hes  (ijl^)  ♦ 16t  (I^*)  (28 

^Green,  A.  E.,  W.  Zema,  Theoretical  Elasticity,  Oxford  at  the 
Clarendon  Press,  1954. 

2Mumaghan,  F.  D.,  Finite  Deformations  of  an  Elastic  Solid,  Dover 
Publications  Inc.,  New  York. 
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Where  X,y 


Lame's  second  order  elastic  constants 


J.,m,n  (Murnaghan's)  third  order  elastic  constants 

q,r,s,t  Fourth  order  elastic  constants 

(Factor  16  is  for  convenience  only) 

Substituting  equation  (27)  into  (28)  we  have 


Vi  = (Ii-3)2 

* [(Il-3)"-2(l2-3)  + 4(Ij-3)] 

+ £ ^ (Il-3)'  + m jj  [(Ii-3)'-3(Ij-3)(l2-3)  + 6(Ii-3)"] 

- ^2-3)] 


+ q(Il-3)'‘  + r[(I^-3)"(l2-3)  - 2(Ij-3)n 

+ s[(Ii-3)(l3-l)  + (Ij-3)"  - (Ij-3)(l2-3)] 

+ t[(l2-3)2  - 4(12-3)  (1^-3)  + 4(Ij-3)*] 

Using  (26),  (13)  and  (5)  and  some  manipulation  we  have 
- 5X1  I •%  ^^1  ^ 5X-i 

'■“p  ■ * “(’b  - * 

+ 2m  Xi(Xi*-l)  ♦ ^ (Xi*-1)*  + 216q(Xi*-l)H3Xi 


(29) 


(30) 


+ 24r(Xj*-l)2(8Xi  - A)  + 8s(Xj*-1)2(2Xi  - i)  ♦ 8t(Xj*-l)*  (TX^  - ^) 
nd 

= T * y(2Xi  CXi*-i)*  ^ (Xi^-i) 

(X  *-l)  + 216q  X (X  *-l)*  ♦ 18r  (X  *-l)*(3X,  - X) 

4Xi  1 ^1  ^ ^ ^ (31) 


+ 2s(Xj*-l)»(Xi  - 
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In  the  above  expression  there  are  two  second  order,  three  third 
order  and  four  fourth  order  elastic  constants.  However,  only  two 
unique  combinations  of  each  third  and  fourth  order  elastic  constants 
appear  in  (30)  and  (31);  these  combinations  can  be  easily  found  by 
collecting  equal  powers  of  the  strain  parameter  These  are 

a = (9Jl+n) 

Third  Order 

B (3il+2m) 

(32) 

Y = (27q  + 9t  + s + 3t) 

Fourth  Order 

6 = (81q  + 24r  + 2s  + 7t) 

Using  (32)  the  expressions  in  (30)  and  (31)  become 

Vp^p  = Xj5(86)  + Xj®(^  + B - 8y  - 166)  + 3y  - | - B + 16y  + 86) 

+ ( - ^ - y + j a - 8y)  (33) 

= Xi*(66-10y)  + Xi'(|  B - 126  + 12y)  + Xj(^  + 2y  - ^ - |b 

1 3X  1 

+ 6y  + 66)  + ( - -2"  - y + 4 a - 8y)  (34) 

The  Equation  of  State 

Using  (7)  and  noting  that  = -Pressure  we  have  the  equation 

-P  = i (3X+2y)  i (Xj*-1)  + a ^ (Xi"-1)*  + ^ (Xj'-l)*  (35) 

The  first  two  terms  on  right  side  of  (35)  represent  Murnaghan's 
equation,  the  additional  term  is  due  to  fourth  order  elastic  constants, 
and  plays  an  important  role  at  high  strains. 
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It  should  be  noted  that  in  Equation  (35)  only  one  combination  of 
each  order  of  elastic  constants  appears.  This  is  the  consequence  of  the 
loniversal  identity  which  has  to  be  satisfied,  independently  of  the 
assumed  elastic  potential.  This  phenomenon  has  far  reaching  consequences, 
as  will  be  explained  in  the  sequel,  in  obtaining  the  equation  of  state, 
purely  from  the  velocity  ratio. 

Numerical  Results  and  Discussion 

As  explained  in  reference  6,  Afp  and  Af^  are  obtained  experimentally 
and  are  related  to  V_  and  V by  Af„  = V_/2d  Af_  = V /2d  where  d is  the 
specimen  thickness.  A spline  function  fit  was  made  to  these  Afp  and 
Af^  and  using  the  universal  identity  (25)  (which  is  independent  of  the 

Q 

elastic  constants  used)  and  Decker's  equation  of  state  , d was  determined 
from  % 


Xi**  dp/dXi  (1*A) 
12po  (Afp*  - 4/3Af3*) 


(1+A  = 1.054)  (36) 


and  Vp,Vg  computed  for  all  pressures. 

£ 9 

The  Frankel,  Rich,  and  Homan  (FRH)  and  Voronov  and  Grigorev  (VG) 
velocity  data  were  both  analyzed.  Since  the  FRH  data  starts  at  25  kbar 
(the  25  kbar  point  having  been  determined  by  matching  velocity  ratios 
of  FRH  to  VG  at  25  kbar)  only  the  VG  velocity  and  3v/8p  were  used  near 
zero  pressure  to  give  the  second  order  elastic  constants  A and  y,  and 
the  third  order  elastic  constants  a and  0.  Having  obtained  the  second 


^Frankel,  J.,  J.  F.  Rich,  C.  G.  Homan,  J.  Geoph.  Res.  Vol.  81,  No.  35, 
6357-6363,  1976. 

®Decker,  D.  L.,  W.  A.  Barrett,  L.  Merrill,  H.  T.  Hall,  and  J.  D.  Barnet, 
"High  Pressure  Calibration:  A Critical  Review,"  J.  Phys.  Chem.,  Ref. 
Date  1,  773-835,  1972. 

^Voronov,  F.  F.,  and  Grigorev,  S.  B.,  "Influence  of  Pressures  up  to  100 
Kbar  on  Elastic  Properties  of  Silver  Sodium  and  Cesium  Chloride," 

Sov.  Phys.  Solid  State,  Vol.  18,  No.  2,  325-328,  1976, 
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and  third  order  elastic  constants  from  VG  alone  the  fourth  order  elastic 


constants  were  obtained  in  two  ways,  by  analyzing  the  FRH  data  and  the 
VG  velocity  data.  In  addition  to  that,  one  of  the  fourth  order  elastic 
constants  y,  was  obtained  by  fitting  the  results  of  the  Decker  equation 
to  our  formulation  of  the  equation  of  state  (Eq.  35).  For  the  ultra- 
sonic data  the  fourth  order  elastic  constants  y and  6 were  computed  by 
minimizing  the  integral  of  the  sum  of  the  squares  of  the  differences 
between  the  experimental  values  of  Vp  and  v^  and  the  theoretical 
expressions  over  the  parameters  y and  6.  The  results  are  given  in 
Table  1. 

In  Figure  3 we  plot  Vp  and  v^  as  a function  of  pressure  for 
inclusion  of  second,  third  and  fourth  order  elastic  constants.  The 
experimental  results  are  also  plotted. 

An  interesting  observation  from  Figure  3 is  the  peak  in  the  shear 
velocity  as  predicted  by  VG  data  using  fourth  order  elastic  constants. 
These  elastic  constants  predict  a shear  velocity  which  peaks  at  about 
80  kbar.  A peak  would  be  interesting  from  lattice  dynamics  consider- 
ations using  the  Bom  criterion  of  phase  changes.^  The  FRH  experimental 
data  does  not  show  such  a peak. 

In  terms  of  our  formulation  the  results  from  the  VG  data  agree 

very  well  with  the  Decker  equation  (see  y in  Table  1 and  Figure  4). 

The  fourth  order  elastic  constants  from  VG  evaluated  up  to  80  kbar 

^Frankel,  J.,  J.  F.  Rich,  C.  G.  Homan,  J.  Geoph.  Res.  Vol.  81,  No.  35, 
6357-6363,  1976. 
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TABLE  1. 


Elastic  constants  obtained  in  our  formulation.  The  second 


and  third  order  constants  were  obtained  from  low  pressure 

9 

velocity  data  of  Voronov  and  Grigorev  . The  fourth  order 
elastic  constants  are  obtained  from  present  analysis  based  on 
velocities  reported  in  ref.  6.  The  fourth  order  constant 
under  source  D was  obtained  from  the  Decker  equation  using 
VG  values  for  second  and  third  order  elastic  constants. 

Source 

VG 
FRH 
D 


ELASTIC  CONSTANTS  (x  lO"  DYNES/CM^) 

Second 

Order 

Third  Order 

Fourth 

Order 

X 

y 

o B 

Y 

5 

1.5430 

1.4699 

-54.8724  -30.845 

4.4701 

18.3613 

II 

11 

II  It 

6.5206 

28.1030 

If 

II 

II 

4.8716 

bprankel,  J.,  J.  F.  Rich,  C.  G.  Homan,  J.  Geoph.  Res.  Vol.  81,  No.  35, 
6357-6363,  1976. 

^Voronov,  F.  F.,  and  Grigorev,  S.  B.,  "Influence  of  Pressures  up  to 
100  Kbar  on  Elastic  Properties  of  Silver  Sodium  and  Cesium  Chloride," 
Sov.  Phys.  Solid  State,  Vol.  18,  No.  2,  325-328,  1976. 


predict  the  equation  of  state  to  270  kbar.  The  equation  of  state 
obtained  from  the  FRH  velocities  is  also  included  in  Figure  4 as  an 
internal  check  since  the  Decker  equation  was  used  to  obtain  the 
velocities.  The  results  from  the  FRH  data  differ  from  the  Decker 
equation  curve  by  about  5 percent  of  the  change  in  the  linear  dimension 
to  270  kbar,  or  on  the  pressure  scale  by  z 10  out  of  270  kbar  (4%). 
Figures  3 and  4 point  up  an  anomaly:  even  though  the  fourth  order 
elastic  constants  fit  the  FRH  deduced  velocities  better  to  270  kbar,  the 
fourth  order  elastic  constants  obtained  from  VG  to  80  kbar  predict  the 
equation  of  state  better. 

Further  work  on  fifth  order  elastic  constants  and  use  of  velocity 
ratios  is  indicated.  Figures  3 and  4 strongly  indicate  that  third  order 
elastic  constants  are  not  sufficient  in  predicting  velocities  or  the 
equation  of  state  beyond  80  kbar.  The  fact  that  the  present  develop- 
ment not  only  derives  the  equation  of  state  to  fourth  order  elastic 
constants  but  enables  us  to  determine  them  from  velocity  information 
will  find  application  in  very  high  pressure  or  high  compression  equation 
of  state  studies. 

In  deriving  the  FRH  results  of  Figures  3 and  4 we  have  used  the 
Decker  equation.  However,  an  independent  equation  of  state  can  also 
be  obtained  if  we  have  good  velocity  data  in  the  neighborhood  of  zero 
pressure. 


i 


i 


I 

I 
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The  computation  is  done  as  follows:  a spline  function  fit  is 
obtained  to  the  available  data  near  zero  pressure  (using  interpolating 
splines)  and  the  derivatives  obtained  at  P = 0 (i.e.  = 1).  From 

(33)  and  (34)  we  then  obtain  X,vi,a,3,  the  second  and  third  order 
elastic  constants.  With  these  we  minimize  the  integral  of  the 
difference  of  the  squares  of  the  experimentally  obtained  Afp/Af^,  with 
respect  to  y and  6 (fourth  order  elastic  constants),  and  the  ratio  of 
equations  (33)  and  (34)  over  the  required  range  of  pressure  (least 
square  fit).  Once  y and  6 are  thus  determined  the  equation  of  state 
is  then  easily  obtained  from  equation  (35).  This  method  has  an 
additional  advantage,  that  it  is  not  necessary  to  measure  the 
deformed  dimensions  of  the  specimen. 

COMPARISON  WITH  AVAILABLE  ELASTIC  CONSTANTS 

The  fourth  order  elastic  constants  are  not  yet  available  in  the 
open  literature.  However,  Barsch^^  has  done  an  extensive  review  and 
analysis  for  third  order  isotropic  elastic  constants.  Comparing  our 
equation  (29)  to  Barsch's  equation  5 we  get  in  his  notation 

n = 7/2  C45g 

“ = ^144  ■ ^456 

I = h (Cj23  2m-n) 


^®Barsch,  G.  R.  "Relations  Between  Third  Order  Elastic  Constants  of  Single 
Crystals  and  Polycrystals,"  Jour.  Appl.  Phys.,  Vol.  39,  8,  3780, 

1968. 
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We  evaluated  data  given  by  Barsch^^  from  four  different  sodium 
chloride  experiments  by  substituting  (37)  for  the  third  order  elastic 
constants  in  (32).  The  results  are  given  in  Table  2.  Since  the 
NaCl  2 data  falls  strongly  out  of  the  pattern  of  the  other  three,  we 
compared  the  a and  & as  given  in  Table  1 with  the  average  of  NaCl  1, 

3 and  4.  The  results  are  in  excellent  agreement. 

Table  1 also  shows  the  satisfactory  agreement  between  the  values 
of  Y and  5 obtained  from  different  experiments,  and  the  Decker 
Equation,  as  discussed  above. 

No  mention  has  been  made  so  far  as  to  the  nature  of  the  elastic 
constants  obtained.  We  note  two  factors  (a)  the  specimen  in  our 
experiment  was  compressed  isothermally;  (b)  velocity  measurements  at 
the  frequencies  of  our  experiment  (~  lOMHZ)  are  adiabatic. 

Hence  in  high  pressure  experiments  we  make  adiabatic  measurements 
on  a state  arrived  at  by  isothermal  deformation.  (The  hydrostatic 
pressurization  is  assumed  thermodynamically  reversible,  which  is  not 
totally  true  for  quasi-hydrostatic  systems.)  The  elastic  constants 
thus  arrived  at  are  known  as  mixed  elastic  constants  in  the  literature. 

No  corrections  have  been  made  from  mixed  to  isothermal  elastic  constants 
because  the  difference  for  NaCl  can  be  considered  small  enough  to  fall 
within  the  uncertainty  of  the  determination.^^ 

^^Barsch,  G.  R.  "Relations  Between  Third  Order  Elastic  Constants  of  Single 
Crystals  and  Polycrystals,"  Jour.  Appl.  Phys.,  Vol.  39,  8,  3780, 

1968. 

llBarsch,  G.  R.  and  Z.  P.  Chang,  "Adiabatic  Isothermal  and  Intermediate 
Pressure  Derivatives  of  the  Elastic  Constants  from  Cubic  Symmetry," 
Phys.  Stat.  Sol.  19,  139  (1967),  and  private  communication  (Chang). 
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TABLE  2.  THIRD  ORDER  CONSTANTS  FROM  FOUR  DIFFERENT  EXPERIMENTS 


AS  GIVEN  BY  BARSCH^®  CONVERTED  TO  OUR  FORMULATION 


Comparison  of  Third  Order  Elastic  Constants  of  NaCl  (x  10^^  dynes/cm^) 


a -56.86  -33.53  -55.76  -51.18 
B -32.33  -24.16  -32.22  -29.96 


Average  of  1,3,4  %Difference  wrt 

Table  1 


-54.27 

-30.76 


1.1 

0.27 


^^Barsch,  G.  R. , "Relations  Between  Third  Order  Elastic  Constants  of 
Single  Crystals  and  Polycrystals,"  Jour.  Appl.  Phys.,  Vol.  39,  8, 
3780,  1968. 
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CONCLUSIONS 


1.  Using  fundamental  principles  we  have  derived  the  equation  of 
state  and  velocities  in  terms  of  second,  third,  and  fourth  order  elastic 
constants,  for  an  isotropic  solid. 

2.  The  results  imply  that  a)  with  sufficiently  accurate  velocity 
data  near  zero  pressure  and  the  velocity  ratio  measured  to  high 
pressures  an  equation  of  state  can  be  determined  outside  of  regions 
where  phase  changes  take  place,  b)  with  sufficiently  good  velocity 
data  near  zero  pressure  or  at  low  pressures  a pressure  scale  cun  be 
predicted  from  experimentally  obtained  velocity  ratio  information. 

3.  The  mechanism  of  computation  can  easily  be  extended  to  fifth 
order  elastic  constants. 

4.  In  general  wave  methods  are  much  more  sensitive  for  the 
measurement  of  elastic  properties  of  matter  than  static  measurements 
because  of  relative  sensitivity  of  velocity  measurements  to  volume 
measurements  at  these  pressures. 

5.  As  a last  point  it  is  useful  to  stress  the  complete  generality 
and  the  high  pressure  predictive  powers  of  the  method  derived  here 
together  with  the  velocity  ratio  concept. 

The  equations  are  derived  from  first  principles  and  apply  to  any 
isotropic  solid  in  a pressure-volume  range  where  no  phase  changes  take 
place.  These  can  be  carried  to  fourth  and  fifth  order  elastic  constants 
and  the  constants  can  be  obtained  to  any  pressure  to  which  the  velocity 
ratio  can  be  measured.  The  equation  of  state  derived  from  these 
principles  with  the  fourth  and  fifth  order  elastic  constants  is 
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expected  to  have  predictive  powers  to  higher  than  the  measured  pressures 
better  than  any  method  heretofore  developed. 

This  will  be  useful  in  obtaining  the  very  high  pressure  molecular 
equation  of  state  for  materials  such  as  rare  gas  solids.  These 
molecular  results  will  then  be  compared  with  available  equations  for 
the  metallic  state  and  the  pressure  of  the  molecular-metallic 
transition  will  then  be  predicted,  together  with  its  estimated 
uncertainty. 
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